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INTRODUCTION 


This  report  describes  the  work  carried  out  durhg 
the  first  three  years  of  research  devoted  to  aerosol 
studies,  both  in  the  UMI3T  laboratories  and  particularly 
at  our  field  station  at  Great  Dun  Fell  (847  m a s 1) 
in  Cumbria.  Of  necessity,  a considerable  amount  of 
time  and  effort  has  been  devoted  to  the  design,  construc- 
tion, installation  and  calibration  of  equipment,  partic- 
ularly at  the  latter  site  and  it  is  considered  worth- 
while describing  the  comprehensive  facilities  which  now 
exist  there  in  some  detail.  This  decription  forms  the 
basis  of  Section  2. 

Some  preliminary  measurements  have  been  made  primarily 
to  ensure  the  accurate  and  reliable  operation  of  the  various 
devices  installed  rather  than  with  any  closely  defined 
specific  scientific  goals  in  mind  at  this  time.  However, 

two  specific  aerosol  projects  have  been  conducted  during 

* 

this  period.  The  former  project  is  a laboratory  investi- 
gation of  the  effects  of  mixing  on  the  evolution  of  wet 
aerosol  spectra,  and  the  latter  project  deals  with  elec- 
trical studies  of  the  turbulent  mixing  of  aerosol  trapped 
at  a temperature  inversion.  These  htudies  are  described 
in  detail  in  sections  3 and  4. 

xhe  final  section  of  this  report  briefly  describes 
and  discusses  recent  work  carried  out  and  likely  future 
developments 
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THE  GREAT  DUN  FELL  RESEARCH  STATION.  AND  ASSOCIATED  FACILITIES 


I 


2.1  IntroducJ ion 

The  UMIST  Field  Research  Station  is  situated  on  the 
summit  of  Great  Dun  Fell,  in  Cumbria,  (Fig  2.1).  The 
facilities  comprise  a stone  building,  situated  847  metres 
above  sea  level  with  adjoining  concrete  forecourt  for  the 
disposition  of  ground-based  instrumentation  and  an  adja- 
cent wooden  mast  upon  which  instruments  can  be  mounted 
up  to  a height  of  20m  above  the  ground.  The  site  is  generally 
accessible  throughout  the  year  by  Land  Rover,  although 
in  extreme  winter  conditions  a snow  cat  is  occasionally 
required.  To  avoid  daily  wastage  of  time  through  travel 
between  the  site  and  UMIST  (a  journey  of  approximately 
2^-  hours)  a caravan  is  located  at  the  foot  of  the  Fell; 
in  addition,  facilities  for  overnight  accommodation  exist 
at  the  station  itself. 

The  building  it  equipped  with  heating,  lighting,  and 
telephone.  A water  tank  and  toilet  have  been  installed 
together  with  cooking  and  other  domestic  facilities.  Thus 
the  station  can  be  continuously  manned  for  periods  ef 
several  days  at  a time. 

A large  laboratory  area  has  been  fitted  out  with 
extensive  bench  space  and  ample  mechanical  and  electronic 
equipment  to  permit  onsite  construction,  modification  or 
maintenance  of  apparatus  when  this  is  necessary.  Travers- 
ing the  laboratory  at  one  end  is  a cloud  duct  which  pro- 
vides a convenient  means  of  sampling  the  external  environmer'4- 


2.i 


A comprehensive  and  flexible  range  of  data  recording 


equipment  has  been  installed  which  we  believe  to  be  capable 
of  meeting  all  our  requirements;  this  facility  includes 
a 12  channel  UV  recorder;  a 15  channel  data  logger;  a 
multichannel  digital  tape  recorder;  three  digideck  recorders; 
and  domestic  cassette  recorders  for  use  in  conjunction 
with  individual  experiments.  The  latter  four  systems  all 
store  data  in  a manner  compatable  with  our  PDP-8e  computer 
at  TJMIST  for  convenient  subsequent  analysis. 

The  Great  Dun  Fell  station,  by  virtue  of  its  location, 
has  seveial  highly  advantageous  features.  Because  of  its 
altitude* it  is  frequently  in  cloud  - for  part  of  approximately 
250  days  each  year.  Also  its  location  and  the  prevailing 
winds  are  such  that  it  is  generally  not  exposed  to  local 
sources  of  pollution.  For  these  reasons  apparatus  has 
been  installed  to  measure  and  record  a wide  range  of  atmos- 
pheric particulates,  both  solid  and  liquid. 

Meusureriont  of  Ait  ken  Nuclei 

The  concentration  of  Aitken  nuclei  is  determined  using 
a photoelectric  nucleus  counter  of  the  type  described  by 
Metnieks  and  Poliak  (1959),  Poliak  and  his  co-workers 
have  carried  out  extensive  investigations  on  the  accuracy 
and  reliability  of  these  devices  and  their  performance 
is  well-known  and  documented.  Two  counters  are  operated 
in  parallel,  Figure  2.2,  and  the  system  is  semi-automated 
by  the  use  of  electrical  solenoid  valves  operated  by  an 
electromechanical  timer.  Air  is  sampled  from  the  cloud  duct, 
either  directly  to  the  counters,  to  determine  the  total 
concentration  of  nuclei  or  via  concentric  cylindrical 
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denuders  across  -which  the  potential  may  be  varied  to 
selectively  descrimirate  against  particles  of  different 
mobilities.  In  this  way  an  indication  of  the  size  distri- 
bution within  the  total  concentration  of  nuclei  may  be 
determined.  Alternatively,  the  denuders  may  be  employed 
with  the  photoelectric  nucleus  counters  to  determine  the 
total  number  of  nuclei  and  the  total  number  of  uncharged 
particles.  Provided  that  the  aerosol  is  in  electrical 
equilibrium  and  that  the  charge  distribution  on  the  popu- 
lation of  aerosol  particles  follows  Boltzmann's  Law  an 
average  size  may  be  determined  for  the  population. 

Electrical  equilibrium  of  the  aerosol  may  be  assured  by 
passing  the  aerosol  over  a weak  radioactive  source. 

The  average  size,  as  deduced  from  the  measured  fraction 
of  the  uncharged  nuclei  in  the  actual  aerosol  is  in  fact 
the  size  of  an  equivalent  monodisperse  aerosol.  (Rich, 

Poliak  and  Metnieks,  1959)* 

The  signals  from  the  counters  are  suitably  processed 
and  recorded  on  the  data  logger,  which  provides  a punched 
tape  output.  The  data  is  returned  to  UMIST  for  subsequent 
analysis  by  a PDP-8e  computer. 

In  principle,  the  photoelectric  nucleus  counter  registers 
all  Aitken  nuclei  (10-  cm  < r < 10-5cm)  upon  which  conden- 
sation may  be  achieved  by  the  supersuturation  resulting 
from  the  expansion  ratio  of  1.21:1.  (A  higher  expansion 
ratio  would  result  in  condensation  on  small  ions.)  However, 
in  practice,  this  size  range  is  limited,  at  the  upper 
end  by  gravitational  settling  of  the  largest  druplets 
formed  in  the  counter,  .and  at  the  lower  end  by  looses  due 
to  diffusion  of  the  smaller  particles  in  the  interconnecting 


1 

tubing  'etween  cloud  duct  and  counter.  Fortunately, 
expressions  for  the  diffusion  of  aerosols  in  laminar 
flow  have  been  found  for  such  situations  (Gormley  and 
Kennedy,  1949),  and  corrections  may  be  made  to  the  raw- 
data  to  allow  for  these  errors.  3y  way  of  illustration, 

_7 

for  particles  of  radius  10  cm  passing  into  the  counter 
at  the  normal  operating  flow  rate  (4  1 min-”'),  the  ratio 
of  particles  reaching  the  counter  to  those  present  at 
the  entrance  n/nQ  « 70$,  and  of  those  particles  reaching 
the  counter  a further  25$  may  be  expected  to  diffuse  to 
the  walls  of  the  counter  prior  to  detection.  Fortunately, 
the  magnitude  of  these  losses  decreases  rapidly  with 
increasing  particle  size,  and  suitable  corrections  to 
the  data  are  conveniently  made  during  data  analysis 
with  the  PDF-8  e . 

Observations  of  aitken  nuclei  concentration  are  most 
conveniently  expressed  as  a cumulative  frequency  distri- 
bution, in  which  the  logarithm  of  particle  concentration 
is  plotted  against  the  percentage  of  observations  less 
than  a given  value  of  particle  concentration.  Representative 
results  are  given  in  figure  2.3.  The  experimental  points 
have  been  omitted  from  this  figure  for  the  sake  of  clarity. 

Preliminary  results  reported  in,  for  example,  the  interim 
report  for  December  1975  have  indicated  the  very  small 
amount  of  scatter  in  the  results.  When  displayed  in  this 
manner  it  is  readily  apparent  that  the  distribution  of 
number  concentration  is  well  approximated  by  a log-normal 
curve  which  provides  a convenient  moans  of  parameterizing 

the  distribution  terms  of  the  geometric  mean,  Z.  and  the 
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These  parameters,  together 
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geometric  standard  deviation  a,  . 

O 

with  the  reinvent  details  of  mean  wind  speed  and  direction 
are  summarised  in  Table  2.1. 
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(curves  A and  J)  ar^  readily  explainable  in  terms  of 
the  incidence  of  anthropogenic  sources  of  pollution. 
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The  only  fixed  local  source  of  pollution  is  a gypsum  mine 
located  8 km  distant  from  the  station  on  a bearing  of 
240°.  Vehicular  activity  or  the  operation  of  the  diesel 
standby  generator  at  GDF  ire  readily  apparent  as  signi- 
ficant departures  from  the  customary  diurnal  variation. 

The  occurence  of  other  short-term  fluctuations  (typically 
an  increase  in  concentration  of  up  to  an  order  of  magni- 
tude for  a period  of  ~ 30  minutes)  has  been  reported  in 
a previous  interim  report.  The  nature  of  the  cause  of 
these  aromoloua  fluctuations  remains  uncertain.  We  intend 


to  continue  the  systemmatic  accumulation  of  data  to  provide 
a statistically  significant  "benchmark"  of  Aitken  nuclei 
concentration . under  various  prevailing  meteorological  conditions 


2.3  Measurement  of  Giant  Dry  Aerosol  by  Inertial  Impaction 


The  measurement  of  number  concentration  and  size 
distribution  of  particles  > 10  cm  is  of  particular  interest 
since  these  particles  are  precisely  those  which  contribute 
to  low  visibilities.  To  investigate  this  size  range  five 
double-stage  impactors  of  the  type  described  by  Jaenicke 
(1971)  were  installed  at  GDF  early  in  the  research  period. 
Much  useful  imformation  on  the  operation  of  the  devices 
was  supplied  by  Dr  R Reiter  of  the  Physikalisch- 
Bioklimatische  Forschungstelle  der  Fraunhof er-Gesellschuf t , 
Garmsich-Partenkirchen,  who  is  currently  working  on  an 
automated  counting  technique  using  an  optoelectronic 

thod.  (counting  Is  performed  at  present  by  conventional 
microscopy).  Ir.  ier  to  t:  • able  to  obtain  accurate  measure- 
ments of  aizt  li.  frit  a*  i ■ i * ir  es  ehti.nl  that  the  col- 


lection effici  r.c.v  ■!'  t ir  lividaal  second  stages  be 


precisely  known.  These  nay  be  calculated  by  consideration 
of  the  aerodynamics  in  the  region  of  the  nozzle,  and  art- 
found  to  be  a sensitive  function  of  J he  mean  particle 
density.  The  variation  of  particle  diameter  at  peak 
collection  efficiency  for  the  five  double  stage  impactors 
for  a range  of  mean  particle  densities  is  given  in  Table  c .2 


Table  2.?:  Calculated  Values  of  Particle  Diameter  of 

Maximum  collection  efficiency  for  various 
mean  particle  densities. 


Mean  Particle 
density  g cm-''. 

Stage  1 

* •_  -*.c 

Particle 
. Stage 

Diameter  d 
2 . Stage  3 

,um 

, Stage  4 

; Stage  5 • 

L « 

, 0.65 

0 . 337 

0.626 

1 .30 

2.79 

, 6.36 

o 

o 

0.233 

0.461 

0.972 

2.09 

! 4.75 

o 

o 

CM 

0.1  46 

0.291 

0.623 

1 .355 

’ 3-05 

With  the  optical  microscopy  techniques  available  to 
us,  we  were  only  able  to  accurately  determine  particle  sizes 
d £ 1 .0  am.  Consequently  by  examination  of  the  deposits 
from  stages  4 and  5 only  were  we  able  to  determine  that 
the  calibration  most  closely  corresponding  to  our  measure  - 
raents  was  for  a mean  particle  density  of  1 .(  g:r.  cm  ' . 

In  order  to  accurately  dote. mine  number  concentration 
it  is  necessary  to  know  the  magnitude  of  the  peak  collec- 
tion efficiency  r.  A comparison  of  both  theoretical  cal- 
culations as  well  as  experimental  determinations  published 
in  the  literature  shows  considerable  discrepancies. 
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Jaenicke  and  Blifford  (1974)  have  reviewed  this  subject, 
and  suggested  likely  causes  for  the  discrepancies.  However, 
in  view  of  the  uncertainties  concerning  the  obtainable 
accuracy  of  number  concentrations  and  size  distributions  of 
atmospheric  aerosols,  together  with  the  inordinately  time- 
consuming  nature  of  the  data  analysis,  it  was  deemed 
advisable  to  postpone  making  these  measurements  until 
the  acquisition  of  the  optical  scattering  particle  counter 
described  in  the  next  section. 

The  theory  of  inertial  deposition  employed  to  deter- 
mine the  collection  efficiency  of  each  sta^e  of  the  iapac- 
tors  demands  that  the  size  of  particle  captured  is  actually 
an  "equivalent  aerodynamic  size",  rather  than  some  speci- 
fic linear  dimension,  whether  real  or  derived,  in  order 
to  account  for  variations  in  parameters  such  as  shape, 
density  and  drag  coefficient  of  the  particles.  It  is, 
therefore,  of  interest  to  examine  the  range  of  shapes  and 
sizes  for  which  a given  impactor  stage  is  descrir inative . 

As  a consequence  of  these  considerations,  it  was  felt 
that  ex  'ruination  of  the  samples  obtained  by  the  imp'  ctors 
by  scanning  electron  microscopy  (HEM)  techniques  could 
be  instructive.  To  this  end,  the  microscope  slides  normally 
employed  were  repl  ced  by  metal  blanks  of  identical  size,, 
at  the  centre  of  which  a hole  was  machined  to  accept  a 
st  :ndard  HEM  sample  stage.  Care  was  taken  to  maintain 
the  face  of  the  HEM  sample  stage  flush  with  the  supporting 
metal  blank  to  minimise  disturbance  of  the  airflow  in 
the  critical  region  of  the  impactor  jet,  so  that  samples 
would  be  represent" tive  of  those  obtained  under  normal 
operating  conditions. 
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Examination  oi'  the  photographs  thus  obtained  (typical 


examples  of  which  were  included  in  the  interim  report 
for  June  1976)  disclosed  considerable  variations  in  particle 
shape  and  size  for  any  given  impactor  stage.  The  valua- 
tions in  overall  size  wore  enhanced  by  the  fact  that  many 
of  the  larger  particles  were  clearly  agglomerates.  Many 
small  particles  at  extreme  m gnii'ication  were  clearly 
crystalline  in  nature.  This  variability  in  the  nature 
of  particles  gives  rise  to  considerable  uncertainty 
in  the  value  of  mean  particle  density. 

2 . 4 Measurement  of  Giant  Aerosol  by  Optical  3c  .ttering 

During  the  latter  part  of  the  research  period  we 
acquired  a ftoyco  Model  225  Aerosol  Particle  Monitor  with 
an  autoscan  digital  display  module  508  and  Model  241  sensor. 
This  device  has  been  installed  it  the  C-DF  field  station 
and  is  currently  undergoing  commissioning  and  calibration. 
The  counter  is  precalibratod  with  monodisperse  sources 
of  latex  particles  in  five  radii  catagories:  0.25-0.35, 
0.35-0.7,  0.7-1 .4,  1. 4-3*0  and  3.0-5.0um.  Since  the  cali- 
bration ie  sensitive  to  the  refractive  index  of  the  aerosol 
being .measured,  we  have  deemed  it  desirable  to  check  the 
calibration  for  natural  aerosol  particles. 

2 . 5 Atmospheric  Pert  iculat...  Coll,  dor 

One  further  device  for  s mpling  dry  atmospheric 
particulates  is  described  briefly  here  for  completeness. 

This  is  the  atmospheric  particulate  collector  designed 
by  A3L  at  WyKR.  In  this  device,  air  to  be  sampled  is 
drawn  up  and  into  the  top  of  the  cylindrical  body 
through  circumferential  orifices  and  then  passes  down 
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the  body  to  a filter  mounted  in  a Gelman  filter  holder 
at  the  ba”e  of  the  body.  The  device  is  aspirated  by  a 
Reciprotor  type  R506  pump  such  that  the  flow  rate  into 
the  intake  is  sufficiently  high  to  overcome  the  sedimenta- 
tion effects  of  gravity  but  not  so  high  that  particle 
inertia  effects  become  serious  enough  to  cause  the  partiales 
to  deviate  from  the  flow  stream. 

Two  such  devices  are  currently  being  operated,  one 
at  UMIST  in  a denc'ly  populated  urban  environment,  and 
the  other  at  GDF. 

.6  Measurement  of  Giant  Wet  Aerosol 

In  order  to  determine  the  size  distribution  of  giant 
wet  aerosol  an  electrostatic-disdrometer  has  been  construc- 
ed  and  tested,  both  in  natural  and  laboratory  conditions. 

The  calibration  programme  has  bee  both  detailec.  and  exten- 
sive in  view  of  the  discrepancies  in  the  behaviour  of 
similar  devices  reported  in  the  scientific  literature. 

A block  diagram  of  the  present  system  is  shown  in 
Figure  2.4.  Essentially,  pulses  from  the  Keily  Probe, 
together  with  the  output  from  the  hot-wire  liquid  water 
content  device,  are  processed  and  recorded  on  a domestic 
cassette  recorder  at  the  field  station.  The  tapes  are 
then  returned  to  UMIST  and  played  back  via  the  n«.  cessary 
ancillary  circuitry,  into  the  PDP-8e  computer.  The  computer 
is  programmed  to  accept  pulses,  sort  into  channels  and 
output  the  total  end  channel  counts  to  t ae  digistore. 

Cloud  liquid  water  content  is  read  at  the  beginnine"  and 
end  of  each  histogram  and  also  output  to  the  digistore. 


The  raw  data  thus  available  from  the  digistore  nay 
then  be  fed  to  the  PDP-8e  which  is  reprogrammed  and  run 
to  output  actual  histograms  of  size  distribution  to  the 


teletype.  The  data  is  available  either  digitally  or  plot- 
ted out  in  suitable  format,  and  provides  in  addition  to 
size  distribution,  total  count,  total  liquid  water  content 
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(ie  T mr  ),  the  mean  radius,  dispersion,  peak  channel  size, 
i 

count  and  liquid  water  content  and  pre-  and  post  run  liquid 
w ter  content  derived  from  the  hot-wire  device.  Typically, 
a single  histogram  is  complied  for  2000  pulses  or  twenty 
seconds  of  data.  Various  system  status  checks  and  moni- 
toring facilities  are  available.  Provision  has  been  made 
in  this  latter  programme  to  incorporate  a correction  factor 
for  each  channel  of  the  histogram  once  the  collection 
efficiency  of  the  keily  Probe  has  been  determined  as  a 
function  of  droplet  size. 

Calibration  of  the  Keily  Probe  has  been  carried  out 
in  those  laboratories  using  a specially  constructed  wind 
tunnel  capable  of  ventilating  the  probe  at  speeds  up  t, 

25  m sec-1 . The  air  is  initially  drawn  through  a humidi- 
fier and  then  droplets  are  introduced  into  the  airstream 
from  a 'spinning  top'  generator. 

The  droplet  spectrum  available  from  this  device  may 
be  selectively  broad  or  narrow.  It  is  determined  by  inser- 
tion of  magnesium  oxide  coated  glass  rods  of  known  collec- 
tion efficiency  into  the  air  stream.  In  this  manner, 
the  preliminary  calibration  curve  presented  in  a previous 
interim  report  has  been  corroborated.  The  integrated 
keily  probe  output  and  the  hot-wire  liquid  water  content 


was  found,  to  bo  excellent;  the  linearity  of  the  output 
from  the  hot-wire  device  was  also  good.  A laboratory 
determination  of  the  collection  efficiency  of  the  probe 
is  currently  being  made. 

2.7  Pollutant  3c-avonginr.  by  Raindrops 

The  distribution,  with  raindrop  size  and  time,  of 
various  pollutants  in  rainfall  is  being  studied  to  gain 
information  about  the  entrainment  of  aerosols  in  rain. 

The  technique  employed  utilises  a collecting  raindrop 
disdrometer  (Georgii  and  Wotzel,  1970)  and  analysis  of 
the  collected  samples  by  electrochemical  techniques  using 
ion-sensitive  electrodes.  Using  this  technique,  concen- 
trations of  selected  ions  as  low  as  1 ppm  may  be  deter- 
mined for  sample  volumes  as  small  as  0.3  ml.  Phis  study 
will  be  augmented  by  precise  determination  of  the  raindrop 
size  distribution  using  a novel  design  of  disdrometer 
developed  at  UKI3T,  specifically  suited  to  obtaining 
accurate  results  in  the  generally  high  prevailing  wind 
conditions  at  G-DF.  In  addition,  electrochemical  analysis 
of  the  samples  obtained  by  the  double  stage  imp  ctors 
is  planned. 
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EVOLUTION  OF  V.'iiT  AdilOSOL  SPIiCTRh 
j . 1 Introduction 

The  evolution  of  droplet  pof  ul  .tior  s formed  by 
condens  tion  in  fo^s  and  clouds  has  been  studied,  experi- 
men  ally  and  theoretically,  by  many  orKers.  The  moot 
detailed  and  reliable  information  on  cumulus  clouds  of 
moderate  depth,  has  been  provided  by  the  airborne  investi- 
gations conducted,  presented  and  dlscusoed  by  arner 
(l%9a,  b,  1970,  1973a,  b). 

In  these  papers  rner  points  to  important  inter- 
rel  ted  featues  of  the  shape  and  development  of  the 
dro  1 t spectrum  - prior  to  the  onset  of  a significant 
contribution  from  coalescence  - which  he  is  unable  to 
explain  quantitatively.  These  are  that  (1)  droplets 
of  the  smallest  deto  table  sine  (diameter  d ~ bum)  are 
found  at  all  levels  within  the  cl-uus,  irrespective  of 
the  distance  from  the  vertical  boindariss;  (2)  the  drop- 
let spectrum  is  consequently  broad,  and  is  commonly 
found  to  be  bimodal;  (3)  the  dispersion  of  the  spec- 
trum - defined  as  the  r-tio  of  the  standard  deviation 


ably  leas  than  1 . liis  papers  provide  convincing  evidence 
for  the  association  of  the  other  spectral  fe  tures  alluded 
to  with  mixing  between  clouu  ana  unaersaturated  environ- 
mental air.  however,  detailed  calculations  ( ' .rner, 1973a) , 
in  which  entrainnmnt  of  undersaturated  environmen  al  air 
was  considered,  yielded  evolving  spectra  quite  unlike 
those  observed.  This  led  arner  to  conclude  that  "simple 
mixing  between  the  cloud  and  the  environment  is  unimport- 
ant is  determining  the  droplet  size  distribution,  at 
least  in  the  early  stages  of  cloud  growth". 

On  the  other  hind,  Mason  and  donas  (1974)  and 
Jonas  and  mason  (1974)  have  presented  culcul  tions,  based 
on  a multi-thermal  model  of  the  air  flow,  from  which  they 
conclude  that  the  major  fe  lures  of  the  droplv  t spectra 
can  be  explained  in  terms  of  mixing  between  the  cloud 
and  its  environment.  This  conclusion  has  been  challenged 
by  V/arner  (1975)  who  argues  that  the  dynamics,  thermo- 
dynamics and  microphysica  ..employed  in  these  calculations 
are  unrealistic.  The  calculati  have  been  defended 
by  Mason  (1975),  but  it  is  clear  that  no  consensus  of 
agreement  exists  as  to  the  reasons  for  the  particular 
features  of  cloud  droplet  spectra  to  which  V/arner  has 
drawn  our  attention. 

In  the  present  work  experiments  are  described  which 
were  designed  to  study  the  effect  on  a population  of 
droplets  formed  by  condensation  of  the  admixture  of 
undone, turated  air. 
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The  appratus  devised  for  the  experiments  is  illus- 
trated in  Figure  (3.1)-  Two  moist  air  streams,  of  temper- 
atures, T^  and  Tg  (TA  > Tg) , were  produced  by  drawing  air 
through  a pair  of  bubblers  partially  filled  with  water. 

A pair  of  centrifugal  traps  was  then  used  to  remove  from 
these  streams  droplets  pr  duced  by  bubble  bursting.  The 
streams  were  subsequently  mixed  at  the  top  of  the  fall 
tube,  and  since  they  were  both  almost  saturated,  a drop- 
let cloud  was  formed  by  condensation.  The  flow  rate  of 
air  through  each  bubbler,  and  Vg,  was  measured  using 
a flowmeter  and  could  be  varied  electrically.  Since  the 
entire  system  was  airtight  the  total  flowrate  through  it 
was  + Vg. 

The  iron  fall  tube  was  cylindrical,  of  length  5m 
and  diameter  0.25m.  Its  interior  walls  were  lined  with 
blotting  paper.  Plastic  honeycombing  was  installed  over 
the  entire  cross-sectional  area  of  the  tube,  near  its 
top  and  its  base,  in  order  to  eliminate  large-scale  turbu- 
lence. This  reduced  the  length  of  tube  over  which  the 
spectral  evolution  could  be  studied  to  about  4m.  The 
minimum  total  flowrate  which  could  be  employed,  about 
80  1 min  1 , corresponded  to  a near,  air  speed  U of  about 
40  mm  s-1  down  the  tube  and  thus  the  development  of  the 
droplet  spectrum  could  be  observed  over  a time  interval 
of  up  to  100  s.  Since  the  largest  droplets  grown  in  the 
tube  were  of  diameter  1 2um  it  is  reasonable  to  assume  that 
their  relative  sedimention  was  negligible  during  their 
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passage  down  the  tube.  Nine  ports,  spaced  at  0.5m 
intervals  within  the  working  section,  were  employed  for 


sampling  the  cloud  droplets,  which  were  collected  on 
glass  slides  coated  with  magnesium  oxide.  The  collec- 
tion efficiencies  of  the  slides  for  droplets  of  the  size 
range  investigated  were  determined  in  a series  of  subsidi- 
ary studies.  Temperatures  wore  measured  at  various  levels 
within  the  fall  tube. 

The  liquid  water  content  within  the  cloud  could  be 
measured  using  an  impactor.  Cloudy  air  was  drawn  into  it 
and  passed  through  a rectangular  nozzle.  Droplets  contained 
in  the  air  were  deposited  on  a slide  covered  with  blotting- 
paper  which  was  subsequently  weighed.  The  values  of  liquid 
water  content,  G,  obtained  in  this  manner  were  self- 
consistent,  and  in  reasonable  agreement , with  those  obtained 
by  integrating  over  the  measured  size  distribution  of  the 
droplets.  The  estimated  uncertainty  in  the  measured  values 
of  C was  + 15 'X'*  It  was  not  possible  to  make  me tsurements 
of  C at  the  same  time  as  sampling  the  droplets  passing 
down  the  fall  tube,  fortunately,  this  deficiency  was 
not  important  since  it  was  found  that  conditions  at  any 
level  within  the  tube  remained  essentially  constant  for 
periods  of  up  to  20  minutes. 

The  effect  on  the  droplet  spectrum  of  mixing  in 
undersaturated  but  humidified  laboratory  air  could 
be  investigated  by  opening  either  one  or  two  circular 
orifices  located  at  or  between  the  level  of  ports  3 or  4, 
(port  1 was  the  highest).  The  relative  humidity  of  this 
air  was  typically  about  90,  • The  diameter  of  the  orifices 
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w ■-! 3 either  6 or  25  mm  and  the  total  flow  rate  through, 
them  (V„)  was  ranged  in  different  experiments,  from  6 to 
40  1 min-'' . In  some  experiments  satur  ted  air  at  a 
temperature  exceeding  that  of  the  descending  cloud  was 
introduced  at  a flow  rate  V~,  through  orifices  located  at 
the  level  of  port  6. 

A number  of  subsidiary  experiments  were  performed 
in  order  to  investigate  the  reproducibility  anu  validity 
of  the  results.  The  homogeneity  of  the  spectrum  across 
the  tube  was  tested  at  various  levels  by  locating  sampl- 
ing slides  at  various  positions  between  the  centre  and 
30  mm  from  a wall.  The  constancy  of  the  spectrum  at  a 
particular  level  and  at  different  positions  at  this  level 
was  studied  by  sampling  at  intervals  during  the  typical 
period  of  an  experiment.  In  all  cases  no  significant 
variations  in  size  distribution  were  found.  Analysis 
of  the  numbers  and  sizes  of  droplets  collected  at  various 
locations  on  the  slides  showed  that  the  distribution  was 
quite  uniform  and  that  the  general  procedure,  which  was 
to  count  and  size  the  droplets  collected  on  a central 
region  of  area  0.6  mm  provided  reliable  spectral  informal  ion . 

A g.od  example  of  the  level  of  reproducibility  attained 
in  these  experiments  is  provided  by  Figure  (3*2)  and  its 
associated  Table  (3.1).  In  this  particular  experiment 
(run  VI)  no  entrainment  of  laboratory  air  occurred  and  tin 
liquid  water  content  C was  measured  at  three  levels  to 
be  0.16  + 0.02  g m 


TABLE  ?«1t  Jaimes  of  mean  drgplot  diameter  cT  and  spectral 

dispersion  a,  for  the  spectra  measured  at  all  nine 
ports  in  Run  VI,  No  mixing, 

Ta=30.5°C;  TB-ia°G;  VA=65  1 min’1;  Vfi=55  1 min"1 
, PORT  *.2  3 4 56789 

I 

| | 

! _ * 

j d (um)  7.2  ^ 7.9  7.5  8.0  8.2  7.9  7.6 

• a j 0.14  0.14  0.13  0.14  0.15  0.13  0.13  0.13  0.1 


The  histograms  and  tables  show  that  although  there  are  con- 
siderable variations#  from  port  to  port,  in  the  number  of 
droplets  located  within  a particular  size  band,  of  width 
1 um,  there  is  reasonable  constancy  in  the  values  #f  the 
overall  spectral  parameters  the  mean,  drop-let-  -diameter, 
and  a,  the  measured  dispersion*  These  parameters  are  defined 
by  the  equations 

I = v N d±/N  (3.1 ) 

i 

a =(£  Ni(di  - d“)2)VN^d  (3.2) 

where  N^  is  the  number  of  droplets  per  unit  volume  in  the 

ith  size  category ? these  have  a diameter  di#  N(=^  N. ) is  the 
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t«tal  drop  concentration.  The  histograms  also  reveal  the 
ab»«ace  pf  systematic  variations  in  the  values  of  N.  This 
finding,  together  with  the  observed  approximate  constancy 
of  water  content  G suggest  that  the  loss  of  droplets  by 
deposition  on  the  walls  wag  negligible.  This  conclusion 
was  confirmed  in  other  experiments.  The  Mean  value  of  N 


for  the  spectra  presented  in  Figure  (3*2)  was  1300  cm  , 

and  the  standard  deviation  was  + 10^.  This  level  of  uncertainty 

was  typical  of  those  obtained  in  the  other  experiments. 

3*3  Results 

In  all  the  acceptable  experimental  runs  in  which 
environmental  air  was  not  introduced  into  the  descending 
cloud  the  mature  spectrum  was  established  by  the  time  it 
reached  port  2,  and  thereafter  remained  essentially  unchanged. 
The  dispersion  a of  the  established  spectra  ranged  from 
0.13  to  0.16,  the  higher  values  being  associated  with 
larger  values  of  N and  C. 

Table  (3*2)  presents  measured  values  of  N,  d,  a 
and  C for  the  experimental  run  VIII,  in  which  undersaturated 
laboratory  air  was  drawn  in  to  the  descending  cloud  at  a 
rate  of  6 1 min-1  through  a hole  of  diameter  6 ram  located 
midway  between  ports  3 and  4.  Table  (3.2)  is  shown  on 
the  following  page.  The  spectrum  existing  at  port  2 is 
very  similar  in  its  properties  to  those  obtained  in  the 
runs  with  no  mixing,  illustrated  in  Figure  (3.2).  The 
spectrum  at  port  3 is  probably  a consequence  of  partial 
dilution  by  the  unders  iturated  stream,  and  is  not  useful. 
However,  the  spectra  at  ports,  4,  5»  6 and  7 are  similar 
in  their  values  of  C and- a> suggesting  that  the  effects 
of  the  mixing  process  are  largly  complete  at  port  4 and 
then  remain  roughly  constant  up  to  port  7.  The  values  of 
the  measured  dispersion  a for  ports  4 and  7 are  mor  ;-or-less 
equal  to  that  at  port  2,  although  the  liquid  water  content 
has  dropped  to  almost  half  its  former  value.  In  the  final 
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Table  3.2:  Values  of  total  droplet  concentration  N,  water  content 


1 

stages  uf  its  journey  through  the  fall-tube  some  growth  of  the 
spectrum  occurs,  presumably  as  a consequence  of  the  observed 
cooling  of  the  descending  cloud.  The  mean  size,  water 
content  and  the  dispersion  of  the  spectrum  all  increase 
between  ports  7 and  9*  If  we  compare  the  spectra  at  ports 
2 and  9 we  see  that  the  overall  effect  of  mixing  in  of  dry 
air  followed  by  some  rcgrowth  of  the  spectrum  is  to  decrease 
C by  1 6^b,  to  decrease  N by  37 $,  to  increase  d by  10$  and 
to  increase  a by  19$. 

It  is  perhaps  appropriate  to  stress  at  this  point 
that  a is  the  measured  dispersion.  Although  condensational 
growth  of  a complete  spectrum  will  always,  in  the  absence 
of  other  factors,  act  to  decrease  the  dispersion,  it  is 
necessary  to  take  account  of  the  fact,  in  our  experiments, 
that  droplets  of  diameter  less  than  6 um  could  not  be 
detected.  The  recorded  increase  in  a and  N when  regrowth 
occurs,  is  a consequence  of  the  growth  to  detectable 
size,  of  droplets  smaller  than  6um  in  diameter. 

The  most  significant  finding,  which  was  confirmed  in 
the  other  mixing  experiments  of  this  type  - irrespective 
of  the  size  of  orifice  through  which  the  undersaturated 
air  was  introduced  - was  that  a is  essentially  unaffected 
by  the  mixing,  even  though  G is  substantially  reduced. 

This  finding  may  be  illustrated  by  comparing  the  spectra 
for  ports  2 and  4 of  run  VIII  displayed  in  the  upper  row 
of  figure  (3*3).  It  is  seen  that  droplets  are  removed 
in  almost  equal  proportions  from  all  size  categories,  and 
the  shape  of  the  spectrum  is  preserved.  There  is  some 
preferential  removal  of  smaller  droplets,  as  evinced 
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by  the  slight  decrease  of  d,  but  the  effect  is  minor. 

Confirmation  of  the  generality  of  this  result  in  these 
experiments  is  provided  by  the  spectra  presented  in  the 
lower  row  of  Figure  (3-3)  and  ports  2,  4 and  7 in  run  XIX. 
In  this  experiment  undersaturated  air  was  introduced  at 
port  3 and  saturated  air,  at  a temperature  in  excess  of 
that  in  the  cloud  at  the  same  level,  at  port  6.  The 
spectral  properties  measured  at  each  port  are  presented 
in  Table  (3.3) • 

Table  3.3:  Values  of  total  droplet  concentration  N,  water  content  C, 
mean  droplet  diameter  d and  dispersion  a for  the  spectra 
me  sured  at  all  sampling  ports  employed  in  the  mixing 
experiment  XIX. 

Ta=39°C;  Tb=19°C;  V =1 20  1 min-1;  VB=150  1 min-1; 


25  1 min-1  ; VB=80  1 min' 
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PORT 
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7 

8 

9 

N | 2300  2550 

(cm-3)  | 

1 500 

1350 

1900 

1800 

1850 

C 

(g  m-3)  | 

► A—  • 

0.64  0.67 

0.35 

0.29 

0.56 

0.56 

0.49 

d ! 

(um) 

7.7  7.7 

7.5 

7.3 

7.9 

8.1 

7,7 

— *- 

nr 

0.17  0.17 

0.16 

0,15 

0.21 

0.20 

0.20 

it 

is  seen,  as  with  run  VIII,  that 

there 

rire  three 

readily 

distinguishable  regimes.  At  port  2,  where  no  mixing  has 

occurred,  the  established  cloud  possesses  a liquid  water 

-3 

content  of  0.b7  g m and  the  measured  dispersion  a is  0.17. 
Between  ports  3 and  6 it  is  found  that  the  admixture  of 
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of  undersaturated  air  has  resulted  in  C falling  to  about 

half  of  its  original  value  with  a and  d being  only  marginally 

changed,  'l'he  saturated  air  introduced  at  port  6 causes 

regrowth  of  the  droplet  population,  and  G rises  to  0.56g  m 

i'he  mean  diameter  d increases  and  the  dispersion  is  • 

increased  to  about  0.20.  i'hus  we  obtain  a similar  r suit 

to  that  of  run  VIII  - the  overall  effect  of  the  mixing 

process  is  to cfecrease  C by  about  1 6$  and  to  increase  a 

by  about  18$. 

Some  simple  calculations  were  performed  in  order  to 
see  how  mixing  might  produce  the  spectral  modifications 
observed  in  run  VIII.  Specifically,  the  objective  was  to 
determine  how  the  spectrum  measured  at  port  2 could  evolve  - 
as  a consequence  of  mixing-  into  that  observed  at  port  4; 
and  also  to  see  if  re growth  of  this  spectrum  coulu  explain 
that  observed  at  port  9*  Various  spectra  and  spectral 
properties  pert  ining  to  these  calculations  are  presented 
in  Table  5-4 
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Table  3.4:  Measured  and  calculated  values  of  observable  droplet 

concentrations  IT  , total  observable  droplet  concentra- 
tion il,  water  content  C,  mean  droplet  diameter  d,  and 
dispersion  a,  for  spectra  at  ports  2,  4,  and  9 in 
the  mixing  experiment  VIII. 

Spectrum  A - measured  at  port  2;  Spectrum  B - 
uniform  dilution  of  A to  take  account  of  mixing 
but  not  evaporation;  Spectrum  G - measured  at  port  4; 
Spectrum  D - calculated  for  Port  4 on  basis  of 
homogeneous  mixing;  Spectrum  E - Calculations  for 
Port  4 on  basis  of  inhomogeneous  mixing  without 
partial  evaporation;  Spectrum  F - Calculated  for 
Port  4 on  basis  of  inhomogeneous  mixing  with  partial 
evaporation;  Spectrum  G - Measured  at  port  9; 

Spectrum  H - C lculated  for  port  9 from  growth  of 
Spectrum  F. 

Note  that  the  droplet  d < 6um  present  in  Spectrum  F are  not 
included  in  the  table,  since  they  could  not  be  observed 
in  the  experiments. 
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19 

0 
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B 
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70 

29 

0 

2 

0 
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0. 

18 

0 
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C 

45 
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47 

22 

6 

1 

0 
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0. 

10 

0. 
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D 
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70 
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0 

2 

0 

0 
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0. 

10 

0. 
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E 

98 
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37 

15 

0 

1 

0 
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10 

0. 
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F 

78 
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47 

19 

0 

1 

0 
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7.5 

0. 

10 

0. 
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G 

42 
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129 

75 

61 

18 

8 

5 
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8.1 

0. 

16 

0. 

178  | 
i 

H 

113 

78 

125 

144 

47 

19 

0 

1 

527 

8.0 

0. 

16 

c. 

178  ! 

As  a consequence  of  introducing  air  of  relative 
humidity  about  90'X  at  a rate  of  6 1 min-1  into  the  des- 
cending cloud,  of  flowrate  122  1 min-1,  the  liquid  water 
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content  is  reduced  from  0*19  g m J at  port  2 (Spectrum  A 
in  Table  (3.4)  to  0.10  g m-^  at  port  4 (Spectrum  C). 

On  the  classical  model  of  mixing,  employed  by  Warner,  Mason 
and  Jonas,  and  others  the  number  of  droplets  per  unit  volume 
will  be  reduced  by  the  fraction  6/122,  which  c 'Uses  N to 
fall  from  840  cm  to  800  cm  - yielding  spectrum  B - 
and  these  remaining  droplets  will  then  evaporate  in  an 
environment  of  uniform  undersaturation.  Although  not 
strictly  correct  , the  form  of  the  evaporation  equation 
used  by  W irner  (1969b)  is  sufficiently  precise  for  our 
purposes,  and  was  therefore  employed.  In  fact  some  simpli- 
fication is  possible  because  the  curvature  and  solubility 
terms  may  be  neglected  to  a first  order  approximation  for 
d > 6 am.  The  evaporation  equation  at  t = 0 may  thus  be 
written 

(r  + a)  dr/dt  = 1 .35  SQ  (3.3) 

where  SQ  (=0.35%)  is  the  undersaturation  which,  when 
exhausted,  will  reduce  C from  0.19  to  0.10  pm.  In  this 
equation,  which  is  appropriate  for  our  measured  temperature 
of  24°C,  r and  characteristic  length  a related  to  the 
condensation  coefficient  are  in  micrometres,  a was  taken 
to  be  5uni,  but  the  precise  choice  is  not  critical.  If 
Spectrum  B is  allowed  to  evaporate,  according  to  the 
equation  (3-3),  with  30  repl  ,ced  by  the  (increasing)  value 
of  S,  until  the  relative  humility  in  the  cloud  has  risen 
to  100%  the  spectrum  D i3  produced.  We  see  that  it  is 
very  different  from  the  Spectrum  C actually  observed  at 
port  4.  The  calculated  value  of  h is  much  larger  and  of 
a and  a much  smaller  thant  the  corresponding  me  isured  vluc.-, 
and  we  conclude  thnt.  the  modifical  ions  to  the  spectra 
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are  inexplicable  in  terms  of  homogeneous  mixing. 

The  crudest  possible  description  of  spectral  modifi- 
cation produced  by  inhomogeneous  mixing  would  be  to  propose 
a process  whereby  the  mixpd  airmass  consists  of  small 
regions  in  which  admixture  of  the  dry  environmental  air 
has  resulted  in  complete  evaporation  of  all  the  droplets, 
interspersed  with  small  jegions  to  which  very  little  ary 
air  has  been  added.  The  average  properties  of  such  a mix- 
ture would  closely  resemble  what  we  observe  ie  a reduction 
of  the  total  count  N and  of  the  liquid  water  content  G but 
no  change  in  the  mean  drop  diameter  d nor  in  the  dispersion 
a.  If  we  assume  that  in  the  regions  to  which  just  a little 
dry  air  has  been  added  some  droplet  evaporation  occurs, 
affecting  primarily  the  smaller  droplets,  wt-  would  expect 
the  average  properties  of  the  mixed  airmass  to  have  a 
slight  reduction  in  N and  increase  in  d compared  to  the 
original  values,  a prediction  which  is  in  close  agreement 
with  our  observed  results. 

As  a rough  illustration  of  this  description  of  the 
role  of  rrixing  on  the  drop  spectrum  - complete  evaporation 
of  some  drops  from  all  size  categori  3,  combined  with  some 
preferential  evaporation  of  the  smaller  ones  - we  note  the 
observations  N = 414  cm  at  port  4 md  N = 52 7 cm  ' in 
the  regrown  spectrum  at  port  9»  Thus  we  assume  th  it  in 
addition  to  the  spectrum  obs< rved  at  port  4 there  were 
113  cm  droplets  below  the  detectable  size  limit.  If,  to 
take  account  of  preferential  evaporation,  we  nrbitarily 
■assume  that  these  unobservable  drops  wer<  of  diar»  ter 
4.0,jm,  and  had  been  removed  in  equal  proportions  from  the 
two  lowest  size  c tagories  (d  --  6,  7um)  we  d riv^  Spectrum  F 
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displayed  in  Table  (3*4).  This  is  seen  to  oe  in  close 
agreement  with  that  observed  (Spectrum  C).  Finally,  if  we 
now  allow  Spectrum  F,  including  the  unobservable  drops, 
to  grow  according  to  equation  (3 ,3)  until  C has  risen  to 
0.16  g m J (and  N to  527  cm  J)  Spectrum  H is  produced.  Its 
values  of  a and  d are  seen  to  agree  well  with  those  of 
Spectrum  G,  observed  at  port  9. 

It  would  be  completely  unjustifiable  to  conclude  that 
these  calculations  - or  the  rudimentary  description  of  tiki 
role  of  mixing  on  which  they  arc  based  - hive  much  quanti- 
tative significance.  However,  it  would  appear  safe  to 
contend  that  the  spectral  modification  observed  ore  quite 
inexplicable  in  terms  of  homogeneous  mixing  but  are  con- 
sistent with  an  inhomogeneous  process  in  \hich  some  drops 
of  all  sizes  are  completely  evaporated,  and  some  overall 
evaporation  of  the  spectrum  also  occurs. 

3.4  Discussion 

As  mentioned  earlier,  W.rner  (1969b)  has  shown  that 
the  dispersion  of  the  droplet  spectra  in  cumulus  clouds  at 
levels  of  up  to  about  200  a is  explicable  when  account  is 
taken  of  the  reduced  rites  of  diffusional  growth  associated 
with  values  of  condensation  coefficient  subst-  nti'  lly  below 
unity.  However,  the  shape  and  dispersion  of  droplet  spectra 
at  higher  levels  requires  the  operation  of  n additional 
process.  Entrainment  of  unaersatur uted  environment nl  air 
is  clearly  the  most  likely  candidate  since  it  produces  both 
the  sub-adiabatic  cloud  water  concentrations  and  - by  acti- 
vation of  fresh  nuclei  - the  continued  presence  of  small 
droplets.  However,  as  previously  discussed,  Warner  (1973a) 
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concludes  from  his  calculations  that  a simple  mixing  pro- 
cess is  unimportant  in  determining  the  drop  size  distri- 
bution. Warner  (1975)  is  also  unable  to  accept  the  multi- 
thermal  model  of  droplet  evolution  presented  by  Mason 
and  Jonas  (1974)  and  Jonas  and  Mason  (1975). 

It  is  recognised  that  the  mixing  experiments  described 
in  Section  (3*3)  were  crude  and  in  essential  respects 
unrepresentative  of  entrainment  within  natural  clouds. 
Nevertheless,  it  is  possible  that  they  focus  attention  upon 
one  important  feature  of  the  mixing  process  - its  inhomo- 
geneity. The  calculations  of  the  cited  workers  have  all 
been  based  on  the  assumption  that  the  mixing  is  uniform  and 
homogeneous.  In  other  words,  they  assume  that 'apart  from 
the  drops  which  are  present  in  the  cloudy  air  replaced  by 
undersaturated  air,  all  droplets  at  a given  level  in  the 
cloud  are,  at  any  time,  exposed  to  identical  conditions 
of  supersaturation  or  undersaturation.  Thus  it  follows 
that  as  they  spectrum  evolves  - and  the  liquid  water  content 
increase  - as  a consequence  of  the  overall  upward  motion 

j 

of  the  cloud  air  the  homogeneous  mixing  process  envisaged 
by  Warner  (1973a)  will  produce  a spectrum  consisting  of 
a gradually  narrowing  peak  at  the  large-diameter  and 
followed  by  a plateau  resulting  from  the  continual 
activation  of  fresh  nuclei.  Such  an  evolution  was  c 
calculated  by  Warner,  and  discarded  by  him  as  being 
inconsistent  with  the  observational  evidence. 
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However,  in  the  pr-  sent  experiments  the  mixing  has 
been  shown  to  be  far  from  homogeneous.  Some  droplets  'ire 
removed,  by  evaporation,  from  11  size  categories  and  - to 
judge  from  the  measured  spectra  and  observed  approximate 
const  incy  of  the  dispersion  - many  other  droplets  are 
negligibly  affected  by  the  mixing  process. 

Although  it  would  be  unjustifiable  to  use  our 
experiment'll  results  to  attempt  a quantitative  treatment 
of  the  effect  of  inhomogeneous  mixing  on  spectral  evolution, 
some  general  predictions  may  be  made.  The  spectra  observed 
n..  :.r  cloud  base  have  been  shown  by  Warner  (1969b)  to  be 
explic  :ble  in  terms  of  condens iticnal  growth  without  the 
intervention  of  mixing.  We  now  consider  inhomogeneous 
mixing  of  this  cloud  air  with  its  environment  as  it  rises. 

V/e  assume  that  the  mixing  occurs  in  the  manner  described 
in  Section  (3*3)  ie  some  droplets  are  removed  from  all  size 
c itogories  and  there  is  also  some  proferenti  1 evapor  ition 
of  the  sm.ll.r  droplets.  In  addition,  following  earlier 
workers,  we  invoke  the  activation  of  fr  sh  nuclei  drawn 
into  the  cloud  by  the  mixing  process.  In  this  situation 
the  original,  large  diameter  peak  will  not  narrow  as  the 
cloudy  air  asoends  since  some  droplets  within  it  will  be 
scarcely  affected  by  the  -axing  process,  while  the  growth 
of  others  - particularly  the  smaller  ones  - will  be 
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significantly  inhibited.  Thus,  as  the  cloud  ascends,  the 
large  diameter  peak  will  broaden  and  gradually  lose  its 
identity.  This  general  pattern  of  evolution  was  commonly 
observed  by  Warner. 

Lnollenberg  (1976)  has  presented  a pair  of  droplet 
spectra  measured  a few  metres  apart  near  cloud  base  which 
appear  to  lend  support  to  the  arguments  presented  herein. 

The  mode  diameter  is  more  or  less  constant  but  one  spectrum 
contains  more  than  twice  as  many  droplets  that  the  other. 

He  states  (private  communication)  that  such  inhomogeneities 
are  often  foud. 

It  should  be  stressed,  however,  that  this  simple  model 
of  inhomogeneous  mixing  has  no  quantitative  support.  Indeed 
it  is  not  susceptible  to  realistic  quantitative  test  until 
more  information  exists  on  tie  nature  of  the  mixing  process  - 
on  all  scales  - within  cumulus  clouds.  For  example,  the 
simple  ideas  presented  in  this  section  do  not  appear  able 
to  explain  the  bimodal  spectra  frequently  (but  not  generally) 
observed  by  Warner. 

An  equivalent  description  of  inhomogeneous  mixing 
can  be  presented  in  terras  of  fluotu- tions  in  sum's  :t:,r.‘i  n 
not  tied  to  fluctuations  in  vertioal  air  velocity. 

Bartlett  (1968)  performed  some  calculations  yielding  the 
conclusion  that  fluctuations  in  supers  duration  resulting 
from  turbulent  variations  in  updraught  velocity  do  not  have 
a crucial  effect  upon  spectral  breadth.  Essentially,  this 
is  because  the  increase  in  sixe  of  a droplet  as  it  rises 
between  two  levels  is  insensitive  to  the  speed  with  which 
it  makes  the  journey  - slower  growth  over  a longer  time 
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in  a lower  supersaturution  (velocity)  air  parcel  just  about 
balancing  faster  growth  over  a shorter  time  in  a higher 


velocity  updraught.  However,  if  inhomogeneities  in  super- 
saturation are  produced  by  the  entr  linraent  of  environ- 
mental air  - as  the  present  experiments  suggest  - spectral 
broadening  can  occur  because  the  velocity/supersatur-  tion 
relationship  is  not  perfectly  reversible. 
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AT  AK  INVERJION 


4 . 1 Ir.tr  oduc  tion 

The  role  of  eddy  diffusion  in  atmosp ^eric  transport 
processes  is  known  to  be  of  gre  it  importance,  but  an 
accurate  quantitative  theory  is  not  available  at  present. 
This  is  largely  a consequence  of  the  difficulty  involved 
in  making  the  required  comprehensive  direct  and  reliable 
measurements,  accordingly,  it  is  appropriate  to  explore 
unconventional  and  seemingly  indirect  techniques  for 
studying  eddy  diffusion. 

We  describe  an  electrical  technique  for  exploring 
eddy  diffusive  mixing  in  the  vicinity  of  a temperature 
inversion.  In  essence,  it  consists  "i  measuring  simul- 
taneously the  vertic-il  electric  field,  E,  and  the  ionic 
conductivity  \ of  the  atmosphere  - which  may  be  subdivided 
into  positive  and  negative  components,  \+  and  Values 

of  eddy  diffusivity  above  and  below  the  inversions,  k1 
and  k ^ respectively,  may  be  determined  by  examination 
of  the  vertical  profiles  of  E and  \ . 

The  feasibility  of  such  an  approach  became  evident 
when  research  flights  in  the  Meteorological  Office  Varsity 
aircraft  revealed  that  the  sharp  decrease  in  \ just  below 
an  inversion  - as  a consequence  of  ionic  immobilisation 
on  trapped  aerosol  particles  - was  not  accompanied  by 
a correspondingly  sharp  increase  in  E.  The  electric 
field  increased  much  more  slowly  with  increasing  distance 
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below  the  inversion  and  its  profile  was  in  good  agreement 
with  that  of  dewpoint.  These  observations,  which  are 
described  in  detail  later  are  indicative  of  the  diffusion 
of  spcae  charge,  originally  created  in  the  region  of  conduc- 
tivity variation,  into  deeper  layers  of  the  atmosphere. 

Much  weaker  evidence  was  found  for  a smaller  but  quali- 
tatively similar  mixing  process  upwards  through  the  inversio 

4 .2  Tr.aar--.-t.ical  Tr  .1::.  ;.t 

he  consider  a temperature  inversion  separating  an 
upper  region  of  the  atmosphere  of  low  aerosol  content 
and  high  electrical  conductivity  A.,  and  a lower  region 
of  high  aerosol  content.  Ions  moving  downwards  through 
the  inversion  become  trapped  and  relatively  immobilised 
on  aerosol  particles  and  the  conductivity  * will  decrease. 
For  simplicity,  we  assume  that  A decreases  linearly  over 
a depth  Z1  until  a lower  oonstant  value  A,,  is  achieved. 

This  assumption  is  well  supported  by  our  experimental 
findings.  If  no  nixing  processes  occur  then  conservation 
of  ionic  air-earth  current  >E  demands  that  E,  has  a 
constant  value,  E^,  above  the  inversion,  increases  through- 
out the  layer  of  decreasing  \,  and  possesses  a constant 
value,  E^  lower  down;  \E(=A . E^ =\ 2&  > ) will  be  constant  at 
all  levels. 

However,  let  us  now  suppose  that  mixing  occurs  between 
these  three  regions,  as  a consequence  of  eddy  diffusion. 
Specifically,  we  consider  the  transport  of  space  charge, 
of  density  p,  created  in  the  region  of  conductivity  , r uii- 
ent,  -md  characterise  this  eddy  diffusion  in  ha- si  regions 
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by  the  diffusion  constants  K.,f  K and  . Charge  trans- 
ported upwards  through  the  inversion  will  constitute  a 
current  moving  in  opposition  to  the  ionic  flow,  while 
that  moving  downwards  through  the  lower  boundary  will 
reinforce  the  ionic  current.  3ir.cc,  in  equilibrium, 
the  total  air-earth  current  J must  be  conserved  at  all 
levels  the  existence  of  this  imugint-a  mixing  process 
will  cause  an  increase  in  XB  around  the  top  of  this  middle 
layer,  and  a decrease  around  the  bottom.  It  is  also 
to  be  expected  that  since  the  vertical  extension  of 
the  field  gradient,  resulting  from  the  transport  of  space- 
charge  from  this  layer,  depends  upon  the  efficiency  of 
the  mixing  process,  it  should  be  possible  to  determine 
values  of  the  diffusion  coefficients  from  the  observed 
structure  of  E.  This  possibility  will  now  be  investigated 


:i  .them. ; ticilly  . 


At  e ch  interf  ice  E and  p must  be  continuous  and 
satisfy  the  equations 

J = \E  + (4.1 ) 

and 


dE  _ -£ 
dz  " eQ 

where  efi  is  the  permittivity  of  free  space, 
equations  may  be  combined  to  give 


(4.2) 
These  two 


J = XE  - e Kd2fi  (4.3) 

° 2 
dz 

The  form  of  solution  of  equation  (4.3)  and  the  meth 
for  establishing  values  of  K1  and  K , will  first  be  illus 
trated  by  considering  a limiting  case  in  which  the  charge 
in  X is  discontinuous  at  z-0  ( ie  z,  = 0,  dx/oz  = . 
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In  these  circumstances  the  solutions  for  h and  o are 


for  z < 


where 


and  for 


o 


where 
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These  expressions  for  £ and  c satisfy  the  equations 
(4.2)  and  (4.3)  and  give  E,  o (and  hence  dB/dz)  continuous 
at  z = 0,  where 


and 


/8  Q /).  4 + / Kj  /j  O 

E = J f — - '■ -■  - 1 

+ /\  ^ K,j 


p = J { 


/\^K^  + / X h 


It  is  seen  from  the  solutions  that  determination  of 
involves  simply  a determination  of  the  exponent  of  n, 
by  curve  fitting  below  the  discontinuity,  and  a knowledge 
of  the  constant  conductivity  \ ^ . An  identical  procedure 
can  be  followed  to  determine  the  eddy  diffusion  coefficient 
above  the  inversion,  K>. 

Having  established  in  principle,  the  feasibility 
of  determining  and  h , from  measurements  of  E and  \ ^ , 
we  now  consider  the  more  realistic  situation,  illustrated 
in  figure  (4.1),  in  which  the  region  of  changing  conduc- 
tivity is  of  finite  depth  7. . 7e  assume  the  conductivity 
gradient  in  this  region  to  be  constant,  and  the  governing 
equation  becomes,  for  C < z < , 


J = p 


§ 0 1 - 
1 1 


X?)1E  - eQK 


d2E 


(4 
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Exact,  solutions  to  equation  (4.4)  have-  been  obtained 

computationally.  However,  a simple  anal;  tic  solution  is 
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available  in  the  case  where  Zy  is  small  compared  with 
the  characteristic  length  describing  the  vertical  dis- 
tance over  which  E changes.  In  this  case  we  may  assume 
that  the  two  terms  contributing  predominantly  to  Jf  which 
must  be  preserved,  are  dx/dz  - which  is  taken  to  be 
constant  - and  eQ Kd"E/dz  . In  this  circumstance  equation 
(4.4)  becomes 


J = Xeo  " eoK7p 
dz 


(4.5) 


where  x = }p  + - XpHz/Zj. 

The  solutions  of  equations  (4.3)  and  (4.5)  are 

T n?z 

z<0,  E = - + Be  ^ (4.6) 

x2 

2 

where  n„  = 


2 eQK2 


0 < z < Z.  , 


E 


X 0z 


or!  - v< 


2 (X1-X?)z? 

+ — ; — ]+Cz,  (4.7) 


z > Z, 


where  n,  = „ 

GoK1 


r ( Z-Z. ) 

E = f + Ae  1 1 

*1 


'1 


(4.8) 


The  constants  A,  B,  C and  EQ  may  be  determined  from  the 


conditions  that  E and  dE/dz  are  continuous  at  z = 0,  Z 


and  we  find 
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This  treatment,  in  common  with  the  exact  solution,  and 
that  presented  earlier  for  the  infinitely  thin  layer, 
predicts  exponentially  varying  fields  E outside  of  the 
layer  of  varying  conductivity.  Thus,  curve  fitting  pro- 
cedures may  be  employed  to  determine  the  exponents  n^ 
and  n0,  and  therefore,  through  equations  (4.8)  and  (4.6), 
values  of  and  • The  application  of  this  procedure 
to  the  analysis  of  the  structures  of  E and  X measured 
during  flights  through  temperature  inversions  is  now 
described  cursorily.  . 

4 . 'jField  ..Studies 

The  field  measurements  were  made  using  the  Meteorological 
Office  Varsity  aircraft,  which  carried:  a cylindrical 
field-mill,  of  the  type  devised  by  Kashmir,  which  measured 
the  vertical  electric  field  E,  and  one  horizontal  comp- 
onent; Gerdian  tubes,  modified  for  airplane  usage,  to 
measure  the  positive  and  negative  tolar  ionic  conductivities, 
\+and  x_  respectively  (x  = x+  + x_);  a Cambridge  hygro- 
meter for  measurement  of  the  dew  or  frost  point,  T^; 
a Hosemount  thermometer  to  record  the  air  temperature,  T: 
and  an  accelerometer.  In  iddition  measurements  were  mad*, 
of  air  speed,  pressure  altitude,  vertical  velocity, 
heading,  ground-speed  and  drift.  All  of  these  parameters 
were  recorded  on  two  3^"  strip-chart  recorders.  Th* 
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electrical  parameters  were  also  recorded  in  digital  fern 
on  magnetic  tape,  in  order  to  facilitate  subsequent  data 
processing-  All  three  recorders  were  synchronised  by 
means  of  a digital  clock. 

Trie  accuracy  of  the  measurements  of  E was  estimated 
to  be  + lOX,  or  5 v m whichever  was  greater,  and  changes 
in  E were  detectable  to  + 1y°.  In  the  case  of  the  ionic 
conductivity  the  absolute  accuracy  of  the  measurements 
was  estimated  to  be  + 15/.  However,  changes  in  conduc- 

—15  _ i _i 

tivity  of  10  n m could  easily  be  detected.  The 
response  times  of  the  field  mill  and  the  conductivity 
tubes  were  made  equal  and  set  at  one  second. 

The  general  flight  procedure  consisted  of  a descent 
to  500'  followed  by  a slow  (500'  per  minute)  ascent  to 
determine  the  heights  of  any  inversions.  This  was 
followed  by  a series  of  "stepped  race-tracks"  such  that 
the  major  inversion  was  passed  through,  either  in  ascent 
or  descent  on  each  of  the  straight  legs  of  the  race-track 
at  rates  varying  from  500  ft  min-1  to  1500  ft  min  '. 

A particularly  clear  cut  and  well  documented  set  of 
measurements  was  provided  by  the  flight  of  14  June  1975- 
The  weather  over  the  UK  was  characterised  by  a surface 
anticyclone  of  1033  nib  centred  over  the  English  Channel. 

The  associated  upper  ridge  was  discernable  up  to  at  least 
100  mb,  and  led  to  a well  marked  inversion  situated  at 
about  880  mb.  The  aircraft  measurements  were  made  over 
the  sea  to  the  south  of  the  Isle  of  Wigh' . 

The  air  temperature  and  dewpoint  recorded  on  the 
initial  climb  from  500  ft  is  shown  in  Figure  (4.J).  Thi. 
top  of  the  aerosol  was  well  rtrkod  visa  illy  at  about 
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3940  ft  above  sea  level*  This  was  confirmed  by  a reduc- 
tion in  turbulence  sensed  by  the  accelerometer  at  this 
height.  The  dewpoint  lapse  rate  suggests  that  a layer 
of  drier  and  presumably  relatively  pollution-free  air 
exist  J at  ibout  700  ft  above  sea  level.  The  sea  surface 
temperature  .t  this  tine  was  approximately  1 3°C  so  large 
seal  ction  is  unlikely  to  have  been  present. 

The  raw  data  points  from  the  field  ana  conductivity 
measurements  avert  ed  over  successive  10  second  intervals 
are  shown  in  Figure  (4.1).  The  variations  in  l and  x 
follow  those  previously  discussed,  the  vertical  depth 
over  which  X cl  mgos  below  the  inversion  being  much  less 
than  that  for  E.  Smoothed  curves,  corrected  for  a zero 
error  in  X,  are  presented  in  Figure  (4.3) > which  also 
displays  t o field-driven  iunic  current  Ex . Although 
EX  is  more-or-less  constant  at  large  distances  from  the 
inversion  it  is  seen  to  increase  and  then  decrease,  mov- 
ing downwards  through  the  discontinuity.  This  indicates 
that  there  is  upward  and  downward  mixing  of  charge  origi- 
nally trapped  on  the  aerosol  particles  in  tie  region  of 
changing  conductivity;  the  downward  mixing  being  much 
more  effective  than  the  upward.  The  theoretical  curves 
based  on  nssu_  ed  diffusion  lengths  of  10  m aoove  and 
100  m below  the  inversion  are  in  good  agreement  with 
those  measurements. 

This  preliminary  study  will  be  extended,  both  theoreti- 
cally and  in  the  field,  where  it  is  planned  to  use  the 
glider  described  in  section  5.3  of  this  report. 
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5 . 1 Electric;!  J in;-  ol~  ;a  rolr 

Experiment  s are  in  preparation  to  investigate  the 
scavenging  of  submicrometre  ;erosol  by  electrified  rain- 
drops. Field  studies  have  shown  that  even  when  the  pie- 
vailing  electric  fields  are  weak  individual  charges  on 
raindrops  may  be  very  substantial.  1‘his  is  because  ther. 
is  often  a mixture  of  positively  and  negatively  charged 
drops,  resulting  possibly  from  their  collision  and  separa- 
tion as  they  fall  to  ground  in  a rainshaft.  Preliminary 
evidence  indicates  that  the  capture  efficiency  is  high, 
even  when  the  aerosol  is  uncharged.  The  experimental 
plan  for  the  first  phase  of  these  studies  is  to  produce 
an  aerosol  of  ki  own  and  narrow  size-range  in  the  band 
0.1  to  5 um  by  evaporating  a population  of  larger  drop- 
lets of  sodium  chloride  solution  produced  by  a spinning- top 
device.  Drops  of  controlled  ch .rge  and  size  will  be 
directed  through  this  aerosol  cloud  and  tin  collection 
efficiency  determined,  for  a wide  range  of  parameters, 
by  measuring  the  amount  of  l.aCl  captured  by  the  drops. 

5 • 2 ,ir,  0;  tic  ;1  :1c  .tt.rin,  ...  a:  . l.a,.  ! ;■ 

A major  instrumental  development  has  been  the  design 
and  construction  of  a novel  device  for  the  measured!,  nt 
of  dropl  t size  distributions  in  cloud; . This  work  was 
c ’nducted  in  two  staffs.  First,  an  instrument  w built 
an  i t jted  in  lal  rat  ry  clouds « Satisf  ctory  j rfor  me 
led  to  the  construction,  for  field  us*  , of  i mor  versa til 
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and  rugged  version,  which  possesses  the  additional  faci- 
lity of  direct  measurement  of  voter  content  and  visibility 
within  clouds. 


The  principle  of  this  technique  is  to  measure  the 
intensity  of  the  light  scattered  at  a set  of  angles  from 
a large  number  of  droplets  contained  in  t he  cloud  - typi- 
cully,  a volume  of  70  cm  is  employed.  A numerical 
inversion  is  then  obtained  of  the  integral  equation  des- 
cribing the  scattering  pattern  as  a superposition  of  the 
individual  patterns  from  all  of  the  droplets  within  the 
scattering  volume.  A matrix  is  used  to  approximate  the 
integral  equation  and  an  inversion  is  obtained  using 
the  method  of  least-squares . 

A requirement  of  a successful  device  based  on  this 
principle  is  that  the  scattering  patterns  are  sensitive 
to  variations  in  droplet  size.  More  precisely,  the 
matrix  of  the  sc  .ttering  functions  should  be  closely 
orthogonal.  Accordingly,  calculations  were  undertaken 
to  determine  the  optimum  set  of  scattering  angles.  These 
led  to  the  choice  of  narrow  angle  forward  scattering, 
which  possesses  the  additional  practical  advantage  of 
convenient  and  accurate  measurement. 

In  essence,  the  instrument  is  an  optical  diffrac- 
tometer with  a horizontal  axis,  and  is  shown  schematic  illy 
in  Figure  5«1.  Illumination  is  provided  by  a three  milli- 
watt helium-neon  laser.  The  scattering  putt  rn  was  measured 
using  a novel  technique.  Fibre  optic  light  guides  were 
embedded  in  a metal  disc  in  order  to  ch  mnel  the  light 
from  the  scattering  pattern  on  to  a photomultiplier. 
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When  the  disc  is  rotated*  a rapid  and  repetitive  scar, 
of  the  scattering  pattern  is  obtained.  The  photomulti- 
plier signal  is  amplified  and  then  digitised,  finally, 
this  numerical  data  is  treated  on  a computer,  and  the 
droplet  size  distributions  are  presented  in  histogram  form. 

We  now  outline  the  extension  of  this  principle  leading 
to  a direct  measurement  of  liquid  water  content.  Calcula- 
tions show  that  the  sum  of  the  intensities  of  the  light 
sc  .ttered  by  a droplet  at  specific  angles  have  been  identi- 
fied which  confirm  this  relationship  to  within  + 2%  far 
droplets  of  r idius  between  1 and  40  am.  This  finding 
enables  the  liquid  water  content  of  the  cloud  to  be  measured 
directly,  continuously  and  accurately.  On  a similar 
principle,  another  combintion  of  scattering  angles  can 
be  selected  such  th  it  the  combined  response  is  proportional 
to  the  surface  area  of  the  droplets.  This  relationship 
is  found  to  hold,  over  the  same  size  range,  to  within 
+ 6>o . Thus  the  reciproc  .1  of  the  measured  signal  may 
be  token  as  an  indicator  of  the  visibility  within  the  clou-. 
These  techniques  for  recording  liquid  water  content  and 
visibility,  in  addition  to  droplet  size  distribution, 
h ve-  been  incorporated  into  the  field-version  of  the 
instrument,  which  has  recnetly  beer,  built. 

The  fibre  optic  device  for  measuring  droplet  size 
distributions  has  been  tested  within  laboratory  clouds. 

Tim  r.  suits  as, reed  well,  over  a wide  ringe  of  conditions, 
with  those  obtained  using  the  st  .ndard  coated  cylinder 
technique s . 
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■■  ri  'n.-nt3  Using  in  Instrumented  ouilplane 

It  is  proposed  to  extend  the  scope  of  the  aerosol 
studies  being  conducted  at  Great  Dun  Fell  b;y  utilising 
an  instrumented  sailpl  me  which  is  capable  of  sustained 
soaring  performance  and  is  also  built  to  withstand  the 
severe  turbulence  which  can  be  encountered  in  large 
convective  clouds.  The  sailplane  (Figure  5.2)  has  a 
wingspan  of  15m,  is  designed  to  fly  wit ain  clouds,  and 
will  be  able  to  carry  about  35  kg  of  instruments.  Its 

_i 

flight  speed  can  range  from  60  to  250  km  hr  1 , its  minimum 
sinking  speed  in  still  air  is  0.7  m sec-1  at  75  tan  hr 
Topics  to  be  investigated  include:  the  measurement 
of  cloud  droplet  spectra  in  cumulus  and  mount  ;in  lee-wave 
clouds,  and  the  vertical  structure  of  the  atmospheric 
aerosol  in  the  lower  3km  of  the  boundary  layer. 

oailplane  measurements  in  the  first  topic  would  comple- 
ment the  existing  ground  b .sed  experiments  at  G-DF  since 
it  would  b<.  possible  to  examine  the  growth  and  subsequent 
modific  ition  of  Iro; 1 t spectra  in  both  cumulus  and 
mountain  wive  cl  v ’s  r.  1 rel  -tt  these  findings  to  other 
parameters  such  is  liquid  water  content  and  height  above 
cloud  base.  The  verticil  velocity  of  the  lir  can  bo 
calcul  i ted  from  the  e-qu  .ti  n of  motion  of  the  sailplane 
(Dye  and  Tootenhoofd,  1973)-  by  centering  the  sailplane 
in  the  strongest,  part  of  the  updraught,  measurements  can 
be  made  whilst  cli.nbin,  towards  the  top  of  the  cloud. 

The  second  topic  is  concerned  with  the  relation  between 
aerosol  concentration  profiles  measured  on  the  sailplane 
and  various  meteorological  parameters  includin<  ’he 
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stability  of  the  layer  and  the  degree  of  turbulent  nixing. 

A particularly  worthwhile  situation  to  examine  is 
thought  to  be  that  of  a layer  of  high  aerosol  concentration 
contained  by  a temperature  inversion  or  hydrolapsc. 

Airborne  measurements  have  already  been  made  of  ion 
concentrations  and  vertical  electrical  field  above  and 
below  temperature  inversions  and  it  has  been  found  that 
the  above-mentioned  electrical  parameters  can  be  used 
to  determine  the  degree  of  eddy  diffusion  occurring 
through  the  inversion.  This  has  been  described  more 
fully  in  section4. 
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